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Abstract. Some rapidly convergent formulae for special values of the Rie- 
mann Zeta function are given. We obtain a generating function formula for 
£(4ra + 3) which generalizes Apery's series for C(3), and appears to give the best 
possible series relations of this type, at least for n < 12. The formula reduces 
to a finite but apparently non-trivial combinatorial identity. The identity is 
equivalent to an interesting new integral evaluation for the central binomial 
coefficient. We outline a new technique for transforming and summing cer- 
tain infinite series. We also derive a beautiful formula which provides strange 
evaluations of a large new class of non-terminating hypergeometric series. It 
should be emphasized that our main results are shown equivalent but are still 
only conjectures. 



1. Introduction 

The Riemann Zeta function is 



oo 1 



k=l 



Apery's formula 1 for ("(3) is 



;i-2) c(3) = -v 



5^(-l) 



fc+i 



Extensive computation has suggested that there is no analogous formula for 
C(5) or CCO- I n other words, if there exist relatively prime integers a and b 
such that 

a ~ (-1)^+! 



f Research supported by NSERC, the Natural Sciences and Engineering Research Coun- 
cil of Canada. 

1 This formula is now commonly associated with Apery's name because it featured in 
his proof of the irrationality of C(3). The formula actually goes back at least as far as [5]. 



2 



BORWEIN & BRADLEY 



then b is astronomically large. Consider however, the following result of 
Koecher [7]: 



fc-i 



<L3) c(5)=2 Sw^gi5<?rgF 

Inspired by this result, the second author, at the suggestion of the first author, 
searched for additional zeta identities of this sort using high-precision arith- 
metic and Maple's lattice-based integer relations algorithms. Within the hour, 
we were rewarded with the following elegant new formula for C(7) : 

Encouraged by this initial success, we searched for and found similar identities 
for C(9), C(11)j C(13), etc. The representation for £(4n + 3) has a convenient 
form in terms of a generating function (1.9), which is our main result (2.1). 
It is curious that there is apparently no analogous generating function in the 
An + 1 case. We refer the reader to the discussion at the end of §8. For now, 
it will be advantageous to exhibit the recursive nature of the formulae in the 
4n + 3 case. 

We denote the power sum symmetric functions P r := Pr S \k) by P := 1 

and 

fc-i 
j=i 

Next, we define a two-place function 

j=l fc=l ^ \k) j=l 

In (1.3) and §9, s = 2 is relevant, but for now we are only interested in the case 
s = 4. Therefore, to minimize symbol clutter we shall occasionally repress the 
superscript, in which case s = 4 should be assumed. With Maple's help, the 
following list was produced: 

|C(3) = A(3,P ), 
^C(7) = A(7,Po) + 5A(3,P!), 
-C(H) = A(11,P ) + 5A(7,P X ) - yA(3,P 2 ) + yA(3,P 1 2 ), 
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-C(15) = A(15, P ) + 5A(11, P 1 ) - - A(7, P 2 ) + y A(7, Pf) 

(1.5) + ^A(3, P 3 ) - ^A(3, P X P 2 ) + ^A(3, J?), 

etc. The first list entry in (1.5) is just a restatement of Apery's formula (1.2) 
and the second list entry is just a restatement of our formula (1.4). From the 
list, it became clear to us that the formula for £(4n + 3) borrows the terms and 
coefficients from the formula for Q(An — 1), except that the first argument of A 
is increased by 4. The number of additional terms is equal to the number of 
partitions of n, and each combination of power sum symmetric functions that 
occurs corresponds to a specific partition of n. Thus, we were led to conjecture 
that 

(1.6) H C (4» + 3) = £ £ E CoF «. (t) , 

j=0 fc=l ^ V k J ahn-j 

where the inner sum is over all partitions a of n — j, the c Q are rational 
numbers to be determined, and if a = (cci, ct2, ■ ■ ■ ) is a partition of n — j, i.e. 
Oil + CK2 + • • • = n — j, then 



pW(fc) := JJpW(*). 



Since it seemed plausible that a generating function could simplify matters, we 
rewrote our conjecture (1.6) in the form 



fc+i 



^yx 4 "c(4n+3) = yx 4 "yy [ ~ \ 2k , y Ca p^\k) 

n=U n=0 J=0 k = l \ k J ohn-i 



j=0 fc = l \ k I ahn—j 

-EE^E^-d) E <^ 4, m 

j'=0 n=j fc = l V re / ahn-j 

OO ^ OO £_-Qjfc + l ^ ^ 

= y xAs y h3 (2k\ i_ x 4 / k 4 y c ° p « 4) w 

s=0 k=l k \k) 1 X /K ahs 



(1-7) =E 



where 

oo 

(1.8) ^(x):=^x s ^c Q P( 4 )(fc). 

s=0 ahs 



Let us denote the number of partitions of a non- negative integer n by p(n). 
By convention, p(0) = 1. We suspected that Ef.{x) had a closed form which 
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might be revealed by determining enough of the coefficients in its power series. 
Fortunately, due to the recursive nature of the formulae we were able to extend 
the list (1.5) without unduly straining Maple's lattice algorithms. This was 
accomplished by introducing only p{n) unknown coefficients for £(4n + 3), 
rather than ^j= i?C?), the actual number of terms involved. Also, when the 
evidence warranted, we supplied the coefficients of as many of the additional 
p(n) terms as we confidently could, based on our ability to recognize patterns 
and extrapolate from previously tabulated values. All identities so obtained 
were subsequently checked numerically - typically to 250 significant digits. 

After having sufficiently extended the list (1.5), we were able to determine 
a good many of the coefficients c a for partitions a of small positive integers, 
and hence the initial terms of the series expansion (1.8). Maple's convert (series, 
ratpoly) feature then produced the following evaluations: 

1, 

1 + 4x 
1 — x 

(l + 4x)(16 + 4x) 
(l-x)(16-x) ' 
(l + 4x)(16 + 4x)(81 + 4x) 
(l-x)(16-x)(81-x) ' 

etc. Thus we were led to conjecture that 

k — 1 .a . 
3 = 1 3 

and hence from (1.7) that 

a..) e **•«*. + 3 »=|e V^t^w n M < i. 

n=0 k=l \k J 1 j=l J 

We restate (1.9) in the next section in the form of a conjectured theorem, and 
discuss some of its implications in the subsequent sections. 

2. A Generating Function Formula for ((4n + 3) 



Ex{x) = 
E 2 {x) = 

E 3 (x) = 
E 4 (x) = 



Theorem A (Conjectured). Let z be a complex number. Then 

" 1 _5f, (-l)^ 1 fi l + ^/j 

£j fc3 (1 _ z 4 /fe 4) 2 ! _ 2 4 /fc 4 1 = 1 i _ zVj 



4/»4 
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Remark. Taking coefficients of z 4 in (2.1) yields our formula (1.4) for C(7)- 
Setting z = in (2.1) yields Apery's formula (1.2) for C(3). In general, taking 
coefficients of z 4n in (2.1) yields a rapidly convergent expansion for £(4n + 3), 
the fcth term of which is a rational function of k whose denominator is a power 
of k times the central binomial coefficient, and whose numerator is a symmetric 
function of partial harmonic sums in 1/j' 4 . 

More precisely, we denote the elementary symmetric functions by 

k— 1 

e { r\h) := lt r } J] {1 + r't) = E • 

J = l l^il<j2<---<Jr<fc-l 

and the complete monomial symmetric functions by 

k— 1 

: = [nn(i-n)" 1 , 

where, as customary, [t r ] means take the coefficient of t r . Then, assuming 
(2.1), we have the following 

Corollary 1 (Equivalent to Conjectured Theorem A). Letn be a positive 
integer. Then 

n oo (_-i\k+l n-j 

(2-2) C(4n + 3) = ^ li) £ 4^iV^)4 4) (*)■ 

7=0 fe=l V fc / r=0 



Proof. Extract the coefficient of z 4n from each side of (2.1). 

In light of the fact that both the complete symmetric functions and the 
elementary symmetric functions can be expressed as rational linear combina- 
tions of the power sum symmetric functions, it is possible to rewrite (2.2) in 
terms of the P a of §1, as in (1.6). However, the formula for the coefficients c a 
appears to be very complicated. Thus, we have replaced the sum over parti- 
tions in (1.6) with a much more manageable sum, at the expense of introducing 
additional symmetric functions into the summand. 

An additional consequence of (2.1) is an attractive formula which provides 
strange evaluations for a large new class of non-terminating hypergeometric 
series. 

Corollary 2 (Equivalent to Conjectured Theorem A). For all positive 
integers n, we have the formula 

,~ „s p / n + 1, n + 1, 2n + in, 2n — in, in, —in 
^ ' ' 6 5 \n+l/2,n,2n + l,n+l + in,n+l-in 
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n 



5 V n / 1 1 4ra 4 + j 4 ' 



Aside. Throughout, we adhere to the standard notation 



F 
p 1 g 



ai, 02, • • • , a p 
bi,b 2 ,...,b q 



E {ai)k{a2)k ■ ■ ■ (a P )k z k 
(bi) k {b2)k---{b q )k kV 



k=0 



where, as customary, 



(a) k := r( " + . fc) =a(a + l)---(a + fc-l). 
Via) 



Proof. We can rewrite (2.1) as a formula for a non-terminating 6 F 5 : 
(2.1) ' : 



1-z 4 

X 6^5 



2, 2, 1 + z + iz, 1 + z — iz, 1 — z + iz, 1 — z — iz 
3/2, 2 + z, 2-z, 2 + iz, 2 - iz 

4 00 



5 S A; 3 (1 - z 4 /fc 4 ) ' 

We note that both sides of (2.4) are meromorphic functions with simple poles 
at z = ±n and z = ±in, where n is a positive integer. We shall see that 
Corollary 2 is a consequence of equating residues of both sides of (2.4) at the 
simple pole z = n. If we denote the requisite residue by R n , then from the 
right side of (2.4), it is clear that 



(2.5) 



Rn — — 



5n 2 



The residue calculation for the left side of (2.4) is more difficult. We have 

fc+i 



Rn 



1 



^ (2) fc (2) fc (1±£± iz) k yr l_ 

k=n — 1 



(3/2), (l) fc (-4) 



"i 3 



4 4 

n — z 
4z 3 



fc+n 



E (2)fc+n-l (2) fc+n _i (1 ± Z ± iz) fc+ra _i -pi- 1^ 
(3/2) fc+n _ 1 (!),+„_! (-4)*-H»-i 1 = 1 j4 _ 



fc=0 
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(2.6) 



(n 4 - z 4 )(2) n _! (2) n _! (1 ± z ± iz) n _! 



4z3( 3 /2) n _! (l) n _! (-4)«-i 



n 



E (n + l) fc (n + l) fc (n ± z zfc zz) fc -j-j- 
(n + l/2) fc (n) fc (-4)fc H 



fc=0 



j=n+l 



(2) n _i (2) n _i (1 ± n ± m) n _i 



n3(3/2) n _! (l) n _i 4 



\ n_1 i 

n)n-l "TT £_ 

n 11 n 4 _ 



X 6^5 



n + 1, n + 1, 2n + m, 2n — m, in, —in 
n + 1/2, n, 2n + 1, n + 1 + in, n + 1 — in 



Comparing (2.5) and (2.6), it follows that 



6-^5 



n + 1, n + 1, In + in, In — in, in, —in 
n + 1/2, n, 2n + 1, n + 1 + in, n + 1 — in 

w4"(3/2) w _i (l) ra _i ^r 1 n 4 -j 4 
5(2) n _ 1 (2) n _ 1 i = i4n 4 +i 4 

2 f4»(l/2) n ^ n(l) B _,g n 4 -j 4 



n! 



(2)„-i f = \4n* + f 



2 {2n\ Vr 1 » 1 ./ 1 



n 



5\n J \^ 4n 4 + j 4 ' 



as required. Thus we have shown that Corollary 2 follows from Conjectured 
Theorem A. That Corollary 2 implies Conjectured Theorem A now follows 
from Mittag-Leffler's Theorem. 

When n = 1, the in Corollary 2 reduces to a 4F3, and we obtain 

Corollary 3. 



4F3 



— ■ —. - / . / 

3/2, 1, 3 



4 
5" 



Aside. Corollary 3 is true, and we have a proof. However, since Corollary 3 
is only a minor consequence of our conjectures, we delay the proof to the end 
of §6, where the proof is used to illustrate some remarks we have to make on 
our methods. 



3. Reduction to a Finite Identity 



As we have said, (2.1) was originally a conjecture based on heavy experimental 
data. However, in the end, we managed to reduce the problem to that of 
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proving a finite combinatorial identity which is beautiful in and of itself, and 
which we have, thus far, been unable to prove. It is 

, , 5^/2£A k 2 S n 4 -j 4 1 

(3 - 1} - 2 ^{k)^Tk^^T¥ = ^ l5SnGZ - 

h — 1 j — 1 

The marvelous connection between the identity (3.1) and the conjectured gen- 
erating function formula (2.1) is presented in the reduction below. 

Reduction. By partial fractions we have for each positive integer k, 

1 tr 1 1 + 4Z 4 /J 4 _ V- Cj(k) 



(3 ' 2) 1 - Z 4 //C 4 II l- Z 4 /? -4 "El 



i - z 4 /fe 4 f = A i - z yf ^ i - z 4 /i 4 ' 



where 



fe-i 



(3.3) c n (k) = ft (1 + 4n 4 /j 4 ) / J] (1 - n 4 //) , 1 < n < k. 

i=i ' j=i 

Substituting (3.2) into the right hand side of (2.1) and interchanging order of 
summation shows that (2.1) is equivalent to 

y 1 _5f 1 y (-l) fc+1 c 3 -(*0 7GC 

2-^3(1-^4) 2^1 Jfe3(») ' 

Clearly, it suffices to prove that for all positive integers n, 

(3.4) E*"(*) = Z3> 



k=n 



where 



(3.5) t n (k) := K J Q , 2fc , l^aeZ. 



5( -l) fc+1 c n (fc) 
2*3 (?) 



Our method of attack is to transform the infinite sum (3.4) into a purely finite 
combinatorial identity. This is accomplished via analytic continuation of the 
summand combined with a process which might aptly be referred to as "Gosper 
reflection" . 

Let n be a fixed positive integer. We wish to extend the definition (3.5) 
to include values of k less than n. One approach is to convert the products 
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implicit in (3.5) into gamma functions. Evidently 
c (k) - lim 1 = xVnA tf 1 + 4 " 4// 

C "W-i™ 1_ x 4 JUL !_ x 4 /(i + 1) 4 

,. r(fc±x±ix)r 4 (fc + l) T 4 (n) r(n + 1 ± x) 
= lim 

^r(i±x±«s) r 4 (fc) r 4 (n + 1) r(fe + i±x) 
r(n + i ± ix) r(i ± g) r(i ± ix) 

T(k + 1 ± ix) r(n ± x) T(n ± ix) 
r(fc±n±in)fc 4 T(2n + 1) 1 r(n + l±in) 

" f(T±n± in) n 4 T(n + k + 1) T(k + 1 - n) T(k + 1 ± in) 
r(n+ 1) (-1)"" 1 r(l ± in) 
X T(2n) (n- 1)! T(n ± in) ' 



(3.6) 



Here, and in the sequel, T(a ± 6 ± c) is shorthand for F(a + 6 + c)T(a + 6 
c)r(a — 6 + c)T(a — 6 — c). It follows that for real k, one can define 



(3.7) 



5e vik T 2 (k + l)T(k±n±in) k T(2n + 1) 
^~~2 r(2fc + l) r(l±n±in)^ 4 r(n + /c + l) 

1 T(n+ 1 ±in) T(n + 1) (-1)"" 1 T(l ±in) 
X r(fc + l-n)r(fc + l±in) T(2n) (n - 1)! T(n ± in) 



Since 1/F(k + 1 — n) = when k is an integer less than n, in view of (3.4) it 
is necessary and sufficient to show that for all positive integers n, 

oo 1 

(3.8) £^ = -3. 

k=o n 

To carry out the reflection process, we need to evaluate t n (k) when A; is a 
negative integer. We shall see that when k is a negative integer, the rather 
forbidding- looking expression for t n (k) given by (3.7) takes a most attractive 
form. From (3.7), 

, , 5 (-1)™" 1 2n(n±in)(±in)n , T 2 (k + 1) 

*n(-l) = -77TZ7Z T\T 7 i I , ; 1 I 1 73 , llm 



2n 4 (n - 1)! (±n ± in)(-l ± n ± in) fc— l r(2fc + l)r(ife + 1 - n) 
5 (-l) n_1 



2n 2 (l + 4n 4 ) (n - 1)! 



x lim r 2 (fc + 2)(2fc + 2)(2fc + l) (fc + 1) ^ 

X ^ (fc + 1)2 r(2A; + 3) f(fc + 2) AJL^ ^ 



n(l + 4n 4 )' 



Let j be a positive integer. One can of course evaluate t n (—j) directly from 
(3.7) by taking the limit as fc — ► —j, just as we evaluated t n (—l) above. 
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However, it is preferrable to introduce the following labour saving device. For 
positive integer k ^ n, define 

m = ~2fc(2fc+l) ((fc+l) 4 -n 4 ) 

an(h t n (k + l) (A; + l) 2 (A; 4 + 4n 4 ) 

For other values of k, define a n (k) by the above expression on the far right- 
hand side. Then for positive integer k, 

t n (-k) = a n (-k)a n (l - k) ■ ■ ■ a n (-2)t n (-l) 

5 -^2j(l-2j)(i-l) 4 -n 4 



n(l + 4n 4 ) 11 (j - l) 2 4n 4 + j A 

j — 2 



n(l + 4n 4 ) 2(k - l)! 2 4ra 4 + k 4 Al 4 n 4 + ; -4 



(3 9) _ 5/2fc\ fc 2 ^n 4 -, 

V ' ; 2n\k ) 4ra 4 + /c 4 ^ 4n 4 + j 4 ' 

Set 

n-l 

p„(A;) := k 2 - j)(/c + j + in)(k + j - in) 

.7 = 1 

<Zn(&) := (A: — n + m)(A; — n — in) 
r n (k) := -2(k + n)(2k-l). 

Then for all integers k, 

t n (k + l) p n (k + l) q n (k) 



4 -4 



t n (k) p n {k) r n (k + l)' 

Furthermore, q n and r n share no linear factors differing by an integer. Ac- 
cording to Gosper's algorithm 2 [4], there exists a polynomial s n of degree no 
greater than 3n — 3 that satisfies 

(3.10) p n (k) = s n (k + l)q n {k) - r n (k)s n (k) 

for all integers k. Define 

r n (k)s n (k)t n (k) 



(3.11) T n (k) := 



Pn(k) 



2 We were led to consider Gosper's algorithm when the first author attempted to get 
Maple to evaluate the sum (3.4) symbolically. Mistyping 'infinity' revealed that Maple could 
evaluate the resulting indefinite sum for specific instances of the parameter n. 
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Using (3.10), it is not hard to show that T n {k + 1) — T n (k) = t n (k) for all 
integers k. Note that since t n (—n) is finite and p n (—n) / = r n (—n), we have 
T n (—n) = 0. It follows that 



m — 1 



T n (m) = ^2 m - 1 ^ -n. 



j = -n 



Also, it is clear from (3.11) and (3.5) that 

lim T n (fc) = 0. 

Thus (3.8) is equivalent to 

- 1 1 
(3.12) T n (0) = £ *«W = -^3- 

j = — n 

Ideally, one would like to prove (3.12) using (3.11). Unfortunately, we do not 
know enough about the polynomials s n to infer the value s n (0) in general. For 
specific values of n, we can use (3.10) to solve for the unknown polynomial s n 
and hence, at least in principle, prove (3.12) for any specific value of n. How- 
ever, using this approach to prove (3.12) in general would require an explicit 
formula for the constant coefficient of the possibly degree 3n — 3 polynomial 
s n . Of course, such a formula can be inferred by assuming (3.12), but to us, 
at least, proving the formula directly seems a formidable task. However, sub- 
stituting (3.9) into (3.12), it is readily apparent that we need only prove the 
beautiful (and apparently non-trivial 3 ) combinatorial identity 

5 ™ (2k\ k 2 ^r 1 n 4 - j 4 1 
( 3 - 13 ) 2^UJ 4n 4 + A: 4 n 4 n 4 +i 4 = ^' " G Z+ 

We discuss the identity (3.13) and some related results in the next section. In 
§6, we examine the process of Gosper reflection in greater detail, where it is 
revealed that identity (3.13) and our conjectured generating function formula 
(2.1) are in fact equivalent. 



3 Professor Zeilberger [10] has kindly informed us that (3.13) and its equivalent hy- 
pergeometric formulation (6.1) do not fall under the scope of identities provable via the 
WZ method. All attempts to prove (3.13) using WZ software, in particular, Zeilberger's 
marvelous package 'EKHAD' and Peter Paule's corresponding Mathematica implementation 
'zb_alg.m', have failed. 
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4. A Combinatorial Identity 



Lemma 4.1 (Equivalent to Conjectured Theorem A). For all positive 
integers n, 



vp-r 5 2 2 /2fc\ 1 TT n4 " 

^ 2 n ^ U J 4n 4 + fc 4 11 4^+ j 4 

fc=i v 7 j=i 

Although we have verified Lemma 4.1 for all positive integers n < 300, we 
have so far been unable to find a proof. The following equivalent proposition 
suggests one possible approach. 

Proposition 4.2 (Equivalent to Conjectured Theorem A). For each 
positive integer n, there exists an even polynomial f n of degree 2n such that 

" x 2 -f 5 (2k\ x 2 k 2 tr 1 x 4 -j 4 

^ x H[ 4iX 4 +j 4 =1 ~ 2^{k J4x* + k* 11 4x 4 +j 4 " 

j=i fe=i v ' j=i 

Clearly, Lemma 4.1 can be obtained from Proposition 4.2 if one sets x = 
±1, ±2, ±3, ■ • • ± n. To see how we arrived at Proposition 4.2, let 

(4.1) ^w:=^v( 2fc )nV+i 2 ) 

V / j=1 

and define a sequence of functions recursively. Put go(x) = 1 for all x and 
for k > let 



fc-i 



r.2 -2 



(4.2, ^w-»w=sniq 

Telescoping (4.2) would prove Lemma 4.1 if we could show that g n (n) = 0. 
Define 



4t 4 

(4.3) fk{x):=gk{x) Yl^ 



4x 4 + j 4 

2~ 



; = 1 x--j 



Then 
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Clearly, g n (n) = if f n {n) is finite. In fact, the evidence strongly suggests 
that each f k is a polynomial. From (4.2) and (4.4) it follows that 

(4.5) (4x 4 + k A )f k _ 1 (x)-(x 2 -k 2 )f k (x) = a k (x), k > 0. 

In particular, (4.2) and (4.5) imply that for all x, fo(x) = 1, fi(x) = 4x 2 — 1, 
f2(x) = 16x 4 + 4, etc. Now Proposition 4.2 is obtained by telescoping (4.2) 
and writing g n in terms of f n . 

We remark that standard telescoping techniques prove the superficially 
similar identity given in 

Proposition 4.3. For each positive integer n, 

5 " fc 4 4 fc ^ n 4 - j 4 

fc = l 7 = 1 J 



Proof. Use 



k— 1 , re— 1 , 



(4.6) J] (a fc _i - b k ) H -L = a - b n [J ^ 



fc=i j=i J j=i J 



with 



4ra 4 + (fc + 1) 4 

a fe = , b k = n - k , fc ^ 0. 

Standard telescoping proves (4.6) for any sequences of a's and b's. In our case, 
we have 

5 

a k -i -b k = -fe 4 , b n = 0, 



and so 



r re fe — 1 4 .4 a d , 1 



4 



4 ti ^4n 4 + (i + l) 

Now cross multiply and obtain 



5v^ kH k k j=l n 4 -/ 



V" - - TT = i 

4 ^ 4n 4 + 1 11 4ra 4 + (j + l) 4 
k=i j=i w y 



from which the claimed identity easily follows. 
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If we try to play the same game using (4.6) to prove Lemma 4.1, it seems 
most natural to define := akin) = 4n 4 + (k + l) 4 for k ^ 0, and then choose 
bk '■= bk(n) so as to satisfy the recursion 



(4 . 7) (at _ 1 _ t0 n|^^(y)n 4n4 -;-; i)4 , *> L 

If we can somehow show that 6 n (n) = 0, then (4.6) implies that 
V -n 2 k 2 ( 2k ) TT re4 ~ J ' 4 = 4n 4 + 1 

2_. 2 n li 4 n 4 + ( J + i) 4 4n +i ' 



A: = l 



which, after cross multiplying, is easily seen to be equivalent to Lemma 4.1. 
Now the recursion (4.7) is equivalent to 



(a k -,{n)- hk{ n))\i^ = \n^( 2k k 

3 = 1 J V 

Thus, b n (n) = is equivalent to 

bj{n) 5 4 /2n\ 

a n -i(n) II —f Ti = -n 

f^- n 4 — j 4 2 \ n / 



3 

i.e. 

n-l 



6o (n) 1 /2n 

LI- 



II „i /; 2 V / 1 J' 

= 1 



which is an equivalent formulation of Lemma 4.1. 

5. An Integral Identity 

Here, we give an exquisite integral evaluation for the central binomial coefficient 
which is equivalent to Lemma 4.1 (3.13) and hence equivalent to our main 
conjecture. 

Corollary 4 (Equivalent to Conjectured Theorem A). For all positive 
integers n, we have the formula 



i r _dv_'jr 



" 1 4y 2 - (j/n) 4 (2n 



j=Q y 2 + U/n) A \n 



Formulae for £(4n + 3) 15 

The proof that Corollary 4 is equivalent to Theorem A relies on the fol- 
lowing conjecture of Wenchang Chu (personal communication): 

Lemma 5.1 (Equivalent to Conjectured Theorem A). For all positive 
integers n, 

n o 2 n ~ 1 / n /o \ 

E^n« 4 +Wn(* 4 -. 4 >=( 2 ;Y 

fc=l j = l ' j = l V 7 



Proof. We'll show that Lemma 5.1 and Lemma 4.1 are inverse pairs. This fact 
is a special case of an inverse pair relationship given in [2] which is equivalent 
to 



(5.1) f[n) = y and n + b n c n Hc k /d k ;n) 

y J , d k Vfc -c fc /4;n + i y ; 

k=r 

if and only if 



where 



fc-i fc-i 
4>(x;k) := JJ(aj + xfy), tp(x;k) := JJ(cj + xd,), 
i=o i=o 



and 



fc-i 



^ m {x;k) := Y[ (cj +xdj 



Setting r = 1, a, = j\ h, = 4, Cj = j\ d 3 = 1, /(fc) = 10fc 2 ( 2 fc fc ) (-l) fc , 
<7(n) = 1/n 2 in the inverse pair (5.1), (5.2) yields the claimed inverse pair 
relationship between Lemma 4.1 and Lemma 5.1. 

We now proceed to show that Corollary 4 is equivalent to Lemma 5.1. 
By a suitable change of variable, the integral identity in Corollary 4 can be 
rewritten in the form 

n— 1 / n /r* 

2n 



(5.3) ^y o ^n(4x 2 -/)/n(- 2 +^ 4 )=(' 



n 
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In view of the partial fraction expansion (3.2), we can rewrite the integrand of 
(5.3), obtaining the equivalent identity 

= (-!) n+1 E^(-)' 

fc=i 

which, in view of the definition (3.3) of the numbers c,t(n), is precisely the 
statement of Lemma 5.1. 



6. Some Remarks on Reflection 



We can rewrite Lemma 4.1 or (3.13) in hypergeometric notation as 

/ 6 -q _p /" 2, 3/2, 1 + n, 1-n, 1 + in, 1 - in _ ^\ 
6 5 y 1, 2 + n + in, 2 + n — in, 2 — n + in, 2 — n — in J 

_ 4n 4 + 1 

5n 2 ' 

an apparently new 4 strange evaluation of a terminating ei 7 ^. We can also 
rewrite (2.1) formula for a non-terminating kR: 



(6.2) 



6-^5 



2, 2, 1 + z + iz, 1 + z — iz, \ — z + iz, 1 — z — iz 
3/2, 2 + z, 2-z, 2 + iz, 2 - iz 



1-z 4 



Observe the dual nature of (6.1) and (6.2). Our process of Gosper reflection has 
taken a non-terminating qF§ at —1/4, and transformed it into a terminating 
6-F 5 at —4, in which certain of the numerator parameters and denominator 
parameters have been exchanged and shifted. 

We can see the dual results of reflection in another way. Let z 4 = — n 4 /4 
in (2.1). The right-hand side terminates, yielding 



(6.3) 



4k 



5A 4 k 
2 2^ ' 



k 



k-1 

n 



rr — j 



^4/c 4 + n 4 2 ^ ( 2k ) n 4 + 4fc 4 11 n 4 + 4j / 
fc=l fc=l V fc / j = l J 

On the other hand, standard techniques show that 
4k 1 



^ 4A; 4 + - 1 



k=l 



n' 1 



2in 2 



il> (1 - in (^±i)) + V (1 + in (^)) - i> (1 + n (±±*)) 



4 See the remarks at the end of this section. 
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-*(¥))} 



1 n 



^ (A; - n/2) 2 + n 2 /4 
Comparing (6.3) and (6.4) yields the following identity: 

5A4^ fc ^ n*-j* _ 1 ^ 1 

2 ( **) n4 + 4fc4 M n4 + 4j ' 4 2n ( fc " ^/ 2 ) 2 + n V4 ' 

Now compare the left-hand sides of (6.5) and (3.13). 

The astute reader will observe a close relationship between the right side of 
(2.3) in Corollary 2 and the summand of Lemma 4.1. In fact, Gosper reflection 
applied to Corollary 2 yields the identity in Lemma 4.1. Since the proof of this 
mirrors the development of §3, we omit the details. We remark however, that 
Gosper reflection easily proves any specific instance of Corollary 2. For the 
sake of brevity, we illustrate this assertion in the case n = 1 i.e. Corollary 3. 
Writing t(k) for the summand of Corollary 3, we have 

= (2) fc (2) fc (±*M-l/4) fc = (fc + l) 2 r(fc±i)r(l/2)(-l/4) fc 
1 ' (3/2)fc(l) fc (l) fc (3) fc r(±i)r(3/2 + k)F(k + 3) ' 

It follows that t(-l) = and t(-k) = for 3 ^ k £ Z. Since 

OO 

E w = °> 

k = — oo 

it follows that 

oo . 

^t(k) = -t(-2) = 

k=0 

which proves Corollary 3. 

To conclude this section, we'd like to offer evidence in support of our 
claim that evaluations (6.1) and Corollary 2 are indeed new. After surveying 
the standard references (eg. [3,1,9]) in the vast hypergeometric literature and 
consulting many of the experts in this area, we have been unable to uncover 
anything remotely like (6.1) or Corollary 2. Hypergeometric summations in 
which the main argument is different from 1 are rare enough. Exceedingly 
rare are summations with complex parameters such as in (6.1) or Corollary 2, 
and neither of our evaluations appears to have a natural generalization. For 
example, there appears to be no generalization of either formula in which in 
is replaced by a general parameter m for example. 
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7. Algorithms and Complexity 



The formulae developed herein lend themselves quite easily to numerical com- 
putation. Algorithms based on Apery's formula (1.2), Koecher's formula (1.3), 
and Bradley's formula (1.4) are particularly simple and are given below. 



a := s + (-l) n+1 /(n 3 c); c := c(4n + 2)/(n + 1); 

od; 

Return(5s/2); 
Zeta5 := proc(d) 

#Compute C(5) to d digits using (1.3). 

N := 1 + [5d/3\ ; Digits := d; a := 0; c := 2; s := 0; 

for n from 1 to N do 

g := 1/n 2 ; s := s + (-l) n+1 (Ag - 5a)/{n 3 c); 

c := c(4n + 2)/(n + 1); a := a + g; 

od; 

Return(s/2); 
Zeta7 := proc(cf) 

#Compute C(7) to d digits using (1.4). 

AT := 1 + [5cZ/3j ; Digits := d; a := 0; b := 0; c := 2; s := 0; 

for n from 1 to N do 

5 := 1/n 4 ; s := s + (-l)" +1 (5a + 5 )/(n 3 c); 

c := c(4n + 2)/(n + 1); a := a + 5; 

od; 

Return(5s/2); 

By Stirling's asymptotic formula for the gamma function, it readily follows 
that 



Zeta3 := proc(cf) 

#Compute £(3) to d digits using (1.2). 

N := 1 + [5d/3\ ; Digits := d; c := 2; s := 0; 



for n from 1 to N do 




fc — > 00, 
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and thus all formulae we have discussed yield 2 binary digits per term asymp- 
totically. Since log 4/ log 10 ~ 3/5, this amounts to slightly better than 1.2 
decimal digits per term. This should be contrasted with the definition (1.1) 
which is asymptotically useless, yielding digits per term. For example, com- 
puting £(3) from the definition (1.1) and applying the integral test to the tail of 
the series shows that the nth tail drops off like 0(l/n 2 ). Thus each successive 
digit requires computing \/T0 times as many terms as its predecessor. To get 
d digits, O(10 d / 2 ) operations are involved. On the other hand, it's not hard 
to see that the algorithms we have presented require only 0(d) operations to 
compute d digits. 

In the preceding discussion, we spoke of operations as operations (such as 
multiplication or division) on numbers. A more realistic evaluation of run times 
must take into account operations on digits. If we take as given that the cost 
of multiplying two d digit numbers is 0(d(logd) log log d), then a crude upper 
bound on the run time for computing d digits using our Apery-like algorithms 
is 

d 

0(^i(log j) log log j) = 0(d 2 (log (i) log log cf). 

i=i 

However, it is possible to adapt these algorithms using the method of Karat- 
suba [6] to yield the highly respectable run time 0(d(\og d) 3 log log d). We coded 
our Apery-like algorithms (sans Karatsuba's optimization) in Maple V (release 
3) and ran them on an Indy R4600PC 100 MHz Silicon Graphics Workstation. 
The following table compares the run times in cpu seconds with Maple's built- 
in implementation of the Riemann zeta function. 

Table of Run Times 





C(3) 


C(5) 


C(7) 


Digits 


Apery-like 


0.4561 


1.8720 


2.8141 


200 


Maple 


8.1720 


8.4600 


8.3462 


200 


Apery-like 


1.1401 


5.5019 


8.0399 


300 


Maple 


28.0742 


28.1819 


28.3860 


300 
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8. Other Dirichlet Series 



For all positive integers n and all real k, let 



577^ 



Then (3.4) becomes 



~ (-l) k + 1 

E pj = !> 1 < n G z 

Thus, for any sequence a\, d2, . . . , we may write 

k=n \ k I 



Let us suppose that n s a„ is absolutely convergent. Summing (8.1) on n 
and interchanging the order of summation, we get 

n=l fc=l V fe / j = l J 

This gives a "formula" for any absolutely convergent Dirichlet series. However, 
(8.2) does not appear to be of much use, except in special cases where we can 
take advantage of known properties of the numbers d n (k). For example, since 

k 

j>(*) = i 

.7 = 1 

for all 1 ^ j ^ k £ Z, putting s = 3 and a\ = a? = ■ ■ ■ = 1 in (8.2) recovers 
Apery's formula (1.2). 

Unfortunately, there seems to be no way to make use of (8.2) or the ideas 
of §3 to obtain a generating function analogue of our result (2.1) for £(4n + 1). 
Since (2.1) started with Apery's formula (1.2) for C(3), one might expect that a 
generating function analogue of (2.1) for £(4n+l) would be based on Koecher's 
formula (1.3) for £(5) and derive from recurrence properties akin to those 
implicit in the list (1.5). However, none of the formulae for £(9) that we 
discovered (and we have good reason to believe there are no others) bears the 
necessary relationship to (1.3). 

We should also point out that even in the An + 3 case, much work re- 
mains to be done, as there are several Apery-like formulae for C(7), C( 11 )j e tc 
which do not arise from our generating function (2.1). In the 4n + 1 case, 



(8.2) 
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the proliferation of formulae appears to be even greater. We have created 
code for systematically listing the formulae for C(13), and ran the code for two 
months or so. The resulting file is over three thousand lines long and contains 
hundreds and hundreds of independent formulae, all having the characteristic 
power of k and central binomial coefficient in the denominator, accompanied 
by harmonic-like sums in the numerator. Classifying the myriad relations and 
interrelations amongst these sums for the various even/odd zeta values would 
be a huge project indeed. 



9. Addendum 



As we later learned, Koecher [8] had given a very simple proof of the following 
generating function for (,{2n + 1), namely 

(9,1) S fc3 (i _ 2/ fc2 ) = ^fe) \2 + l-z 2 /k 2 

If n is a non-negative integer, extracting the coefficient of z n from each side of 

(9.1) produces the formula 

(9.2) C(2n + 3) = ^f;^-^(-ire( 2 )(fc) 

k=l \ k I 

n 00 / i\k+l 

+2EES»(-irei ! ;- ( *), 

j — 1 k — 1 V k ) 

from which (1.2) and (1.3) follow as special cases. Here, the el s \k) are the 
elementary symmetric functions defined in §2. 

Despite the fact that Koecher 's generating function (9.1) gives formulae 
for all odd Zeta values, there is a very real sense in which (9.1) is inferior to 
our generating function (2.1). In (9.1), among other things, the fourth powers 
that feature in (2.1) are replaced by squares. This results in redundant terms 
in his zeta formula (9.2) for n > 1. For example, n = 2 in (9.2) yields 

( 9 - 3 ) C(7) = 2gL^ 



which should be compared with our more compact formula (1.4). To enable a 
more detailed comparison, we rewrite (9.3) in the notation of §1. Then (9.3) 



) n (! - - 2 /i 2 ) ■ 



k-l 



k=l K U ) , = 1 J 



2 ^ kH 2 , k ) ^ Pi 2 ' 

k = l ^ \k I l<j<l<k-l J 
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becomes 

(9.4) C(7) = 2A(7,P (2) ) - 2A(5,P 1 (2) ) + |a(3, p[ 2) p[ 2) ) - {^A(3,P 2 (2) ), 
whereas (1.4) is simply 

(9.5) C(7) = §A(7,P W) + f A(3,PW). 

Since Pq 4 ^ = Pq 2 ^ = 1 and P^ = P^ , we see that the middle two terms of 
(9.4) are redundant. Indeed, lattice-based reduction shows that 

2A(7,P (2) ) + 8A(5,P 1 (2) ) - 5A(3,P 1 (2) P 1 (2) ) + 55A(3,P 2 (2) ) = 0. 

As far as we can tell, in contrast with the formulae derived from (9.2), there 
are no redundant terms in our formulae for £(4w + 3) which come from (2.1), at 
least for n < 12. It goes without saying that, despite our best efforts, Koecher's 
proof of (9.1) apparently cannot be adapted to prove (2.1). It seems that (1.4), 
and more generally (2.1), is a much deeper result. We should also point out 
that merely bisecting Koecher's generating function (9.1) will not yield (2.1), 
nor any new zeta formulae. 

References 

[1] W. N. Bailey, Generalized Hypergeometric Series, Stechert-Hafner, New 
York, 1964. 

[2] Wcnchang Chu, Inversion Techniques and Combinatorial Identities, Bol- 

lettino U.M.I. 7-B (1993), 737-760. 
[3] Ira Gessel and Dennis Stanton, Strange Evaluations of Hypergeometric 

Series, Siam J. Math. Anal., Vol. 13, No. 2 (March 1982), 295-308. 
[4] Ronald L. Graham, Donald E. Knuth, and Oren Patashnik, Concrete 

Mathematics, Addison-Wesley, New York, 1989, pp. 224-228. 
[5] Margrethe Munthe Hjortnaes, Overf<pring av rekken ^fc°=i(V^ 3 ) til et 

bestemt integral, Proc. 12th Cong. Scand. Maths, Lund 10-15 (Aug. 

1953 (Lund 1954)). 

[6] E. A. Karatsuba, Fast Evaluation of C(3), Problemy Peredachi Informat- 

sii, Vol. 29, No.l (1993), 68-73. 
[7] Max Koecher, Klassische Elementaire Analysis, Birkhauser, Boston, 1987. 
[8] Max Koecher, Letter, Math. Intelligencer 2 (1980), 62-64. 
[9] L. J. Slater, Generalized Hypergeometric Functions, Cambridge Univ. 

Press, London, 1966. 
[10] Doron Zeilberger, Personal e-mail, May 10, 1996. 



Formulae for £(4n + 3) 



Jonathan Borwein &; 

David Bradley, 

Centre for Experimental 

and Constructive Mathematics, 

Simon Fraser University, 

Burnaby B.C., 

Canada V5A 1S6. 

jborwein@cecm.sfu.ca 

dbradley@cecm.sfu.ca 



